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ABSTRACT 

In this paper a theoretical model is developed to 
allow for the calculation of the synchrotron emission 
arising from high energy electrons trapped in the Van 
Allen belts of a planet with a dipole magnetic field. 

The model is general enough to allow for the calculation 
of the intensity of radiation received by an observer at 
any distance from and any latitude about the planet. 

The model is used to compute the emission from the 
Earth’s Van Allen belts that one should expect at various 
latitudes at a distance of 1.92 Earth radii, the position 
of the Radio Astronomy Explorer (RAE) satellite that was 
launched in 1968 [12], for the frequencies 1.3 MHz and 
2.2 MHz. The distribution of high energy electrons in the 
Earth’s Van Allen belts used in the model was obtained from 
data published by Vette, Lucero and Wright [ll]. 

The emission predicted by the theoretical model is 
compared with data on the distribution of radiation from 
the Earth’s Van Allen belts gathered by the RAE satellite. 
These comparisons show that emissions from other regions 

of the magnetosphere interfere with the observance of the belts, 
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I . INTRODUCTION 


Since the discovery, by Van Allen [8] of the exist ance 
of the region above the Earth where there exist charged 
particles trapped in the magnetosphere, several attempts 
have been made to develope models that would allow -one 
to compute the electromagnetic radiation emitted by these 
trapped charged particles (electrons) . 

These earlier models (by Vesecky [10] for example) 
would only allow one to compute the emission from the Van 
Allen belts in the geomagnetic equator. Also, when the 
earlier models were developed there was no available data 
on observations of the emission to enable one to make 
comparisons of the observations with the models’ predicted 
values of the emission from the belts. 

The model developed in this paper permits the computation 
of the emission from the belts one should observe at any 
position about the Earth. Comparisons of the model's 
predictions are also made with the data on the emissions 
from the belts gathered by the RAE (Radio Astronomy Explorer) 
satellite launched in July of 1968 [13]. 

II. THE EARTH’S MAGNETIC FIELD 

The Earth's magnetic field can be approximated by a 
dipole field out to several Earth radii. Since we shall only 
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be interested in the lower magnetosphere, (out to approxima- 
tely two Earth radii) where there are high energy trapped 
electrons, the dipole approximation of the Earth's magnetic 
field shall be used in our model. 

The strength of a dipole magnetic field is given in 
the equation below as 

M 2 1/2 

B(r,X e ) = ~ [1. + 3. sin X g ] (1) 

r 

where r is the distance from the center of the dipole, X g 
is the magnetic latitude and M is the magnetic moment. 

III. THE VAN ALLEN RADIATION BELTS 

We now have a pretty complete picture of the radiation 
belts. We know now that the energetic particle distribution 
is made up of high energy electrons in the inner zone and 
energetic electrons in the outer zone. Lower energy particles 
have also been mapped. Particles in the outer zone especially, 
the electrons, show large time variations. The electron 
flux can vary more than an order of magnitude in less than 
a day, responding strongly to magnetic disturbances. 

IV. THE DEVELOPMENT OF THE MODEL 

a. The Coordinate System 

The system of coordinates used in this paper are 
geomagnetic coordinates. The origin of this coordinate system 
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is at the center of the magnetic dipole. For the co- 
ordinates of the observer, we shall use r ,X q and 0 q 

where r is the radial distance to the observer, X is 
o o 

the latitude of the observer and 0 is the longitude of 

o 

the observer. For the coordinates of the radiating 

electron we shall use r ,X and 0 where r is the radial 

e e e e 

distance to the radiating electron, X is the latitude 

© 

of the radiating electron and 0 is its longitude. Since 

the emitting electron is confined to move along a particular 

field line it is more convenient to express r in terms 

e 

of that field line. For a dipole field the equation of 
a field line is 

r = re cos^ X (2) 

e o e 

where re Q is the distance from the center of the dipole to 

the point where the field line crosses the dipole equator. 

However, re Q can be expressed in terms of planetary radii, 

R where re = LR . Thus, L called the L-shell can be used 
E o E 

to define the field line of the gyrating electron and in 

most places in this paper we shall use the coordinates L, 

X ,0 instead of r ,X ,0 . 
e e e e e 
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b. The Emission from a Gyrating Electron in a 
Magnetic Field 

The differential rate at which energy is emitted 

spontaneously per unit solid angle per frequency interval 

du> , called the coefficient of spontaneous emission T) , by 

an electron gyrating about a magnetic field line is [l]. 

2 
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where e is the electronic charge, e Q is the permittivity 

of free space, C is the speed of light, u) is the observing 

frequency, I^Cx) is a Bessel function of order m and argu- 
/ 

ment X, J (x) = dJ (x)/dx and 9 is the angle between the 
mm 

direction of the emitted radiation and the magnetic field 
direction. The quantities Py and P^ are 

P II = 0 cos a and P^ = P sin a (4) 

where P = V/C and a is the pitch angle of the gyrating 
electron. It is the angle between the particle's velocity 

— f 

vector V and the direction of the magnetic field. The 
quantity m is also the harmonic number for which radiation 
is being emitted, m = 1,2,3... The argument of the Bessel 
function, x is 

X = u>/u> P sin 8 (5) 

O X 
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and the argument of the delta function, Y is 

Y = mu) - » (1 - B n cos 8) (6) 

the quantity u> is defined as 

% = \ 1-02 (7) 

where ur , the gyrofrequency , is 

u) b = eB/m (8) 

Assuming that the emission is not affected by the in- 
tevening medium between the observer and the radiating 
electron, the angle 0, in terms of geomagnetic coordinates 
can be found from the expression 

cos 8 = sin y sin X - cos y cos 9 e cos X cos 0 (9) 

- cos y sin 0 Q cos X sin 0 
© 


where 
cos 0 

sin 8 


sin X 


cos X 


r cos X - t cos X cos 8 
o o e e e 


, 2 2 2 2 „ 1/2 
(r cos x + r cos X - 2r r cosx cosx cos 8 ) 
e eo o oeoee 


r cos X sin 0 
e e e 


. 2 2 22 v 1/2 
(r cos X + r cos X - 2r r cosX cosX cos0 ) 
e eo o oeoee 


r sin X - r sin X 
o o e e 


2 2 1/2 
(r + r - 2r r (cosx cosX cos 8 + sinX sinX )) 

o e oe oee oe 


r cos X - r cos X 


2 2 „ . „ 1/2 
(r + r -2r r (cosx cosX cos8 + sinX sinX )) 
o e oe oee o e 
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sin Y 


2 2 
cos X - 2. sin X 0 
e £ 

2 . 1/2 
(1+3 sin X ) 

© 


and 

cos X sin X 0 

6 a 

COS Y = o v/o 

(1 + 3 sin 2 X 2 ) 7 

c. The Particle Distribution Function 

For one electron, if we assume that no radiation 

is absorbed by the medium, the power per unit area for a 

particular frequency or flux P^ is 

r\ (X ,L ,a , E, 0) 

The value of |r| is, in terms of geomagnetic co- 


( 10 ) 


ordinates; 


1/2 


rl= r +r - 2r r (cosx cosX cos0 + sinX sinX 
1 o e oe o e e o e 


( 11 ) 


For more than one electron, to obtain the total flux 
received by an observer, we must sum the individual contribu- 
tions of each particle or electron. Let us assume that the 
radiating medium is continuous. This allows us to write 


dP, = P w (X e ,L,r,a,E,9) dN 


( 12 ) 


where N is the total number of particles in the medium. 


But we can write dN as, 

dN = p (X ,9 ,L,a,E) dE da SV (13) 

e e 

where p is the electron density distribution function and 
dV is a differential volume element of the medium. In terms 


of geomagnetic coordinates 
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( 14 ) 


2 3 7 

dV = L cos X dX d0 dL 

E e e e 

where R is the radius of the Earth, L is the L-shell 
E 

and X and 9 are the geomagnetic latitude and longitude 
e e 

respectively. 

We now have for the differential power 

dP = P (E,X ,9 ,L ,a , r ) p (X ,9 ,L,a,E) X 
u) a) * e e o e e 

XL 2 R 3 cos 7 X dX d9 dL dE da 
E e e e 


(15) 


Considering the particle distribution function we shall 
assume azmutal symmetry by eleminating 9 e from the distribu- 
tion function. This gives 

p = p (X ,L ,a , E) (16) 

© 

According to D. Chang [2], for a dipole magnetic field, 
considering Liouville's theorem 

2 .1/4 

(1 + 3 sin X ) 

p (X ,L ,a , E) = 3 p(0°,L,E,a E ) (17) 

cos X 

e 

where p(0°,L,E,a ) is the electron distribution at the dipole 

E 

equator. The variable a £ is the pitch angle of the electron 
at the equator. 

Let us write p (0°,L,E,a ) _ as 

p( 0 °,L,E,a E ) = p(L,E) f^) (18) 

and choose a value of f(a ) in accordance with D. Chang [2] 

l!i 

in his work on radiation from the planet Jupiter, that is 
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(19) 


f( V " Sin a E 


where P is any positive real number. Thus for our distribu- 


tion function we have 

(1. + 3. sin 2 X ) 1/4 

p(\ e ,L,a E ,E) = — sin P a p(L,E) 

cos X E 

e 

Now, since sin a = sina cos 3 X / (1 . + 3. sin 2 x ) 1/4 we 

e e ’ 


have; 


P (X ,L,a,E) = 
© 


(cos X ) 
e 


3 (P-1) 


(1+3 sin 2 X ) 1/4(p-1 > 
e 


sin a p (L, E) 


( 20 ) 


( 21 ) 


For the limits on the pitch angle distribution we 
want and a , the limits on a , to be such that for a 
given X ^ and L,^ will be the minimum value of a that will 
allow the electron to mirror before crashing into the Earth 
and a 2 to be the maximum value of a that will allow the 
electron to mirror before crashing into the Earth. 

For two point a and b on a dipole field line 


sin a (1+3 sin 2 x )"^ 2 cos^x 

a _ ea eb 

. 2 2 1/2 6 
sm a (1+3 sin X ) cos X 

b eb ea 

If we let the particle mirror at X the latitude where 

eb 

the field line enters the surface of the Earth, then, we 
must have that 

Vo 1 

sin = sin z = 1 and cos X , = — 

b eb a : 

Letting a = a and X = X , we get that 
m a e ea 


( 22 ) 


( 23 ) 
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r „ 2 .1/4 

[1. + 3. sin X ] 

sin a = o — (24 

m L /2 [4. - 3./L] cos 3 X e 

It can easily be shown that a_ = n-a , since the particle 

Z m 

with pitch angle on its return trip after mirroring, 

will have changed directions by 180 degrees. 

Letting P = 2 in equation (21) will give the best 

description of the pitch angle distribution of electrons 

along a magnetic field line in .the Earth's lower magnetosphere 

as experimental evidence has shown [4], Thus we shall use 
2 

f(a) = sin a throughout the paper. 

We finally have for the distribution of electrons 

3 

cos X 2 

P (X ,L,a , E) = r TJ 2 sin a p(L,E) (2i 

e (1+3 sin X ) 

e 

where p(L,E) will be determined by using experimental data. 


d. The Limits of Observation of the Radiating Region 
of a Planet 


Before proceeding to make computations of radiation 
from trapped electrons in a planet's dipole magnetic field, 
we must consider the fact that part of the radiating medium 
will be obscured by the planet. This being the case, we 
must allow for this fact when we make integrations over 


the entire observable radiating medium. The limits of 
integration over the latitude X and longitude 8 are a 
function of the position of the observer with respect to 


the planet. 
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In order for an observer to receive the emitted 
radiation from a particular section of the radiating 
medium about a planet, the condition 


cos X cos X cos (0 - 0 ) + sin X sin X^ 

o e o e o e 


r r 
o e 


(r 2 - R?) (r 2 - R 2 ) 
( o E e E 


1/2 


must be satisfied, where r , X and 9 are the position 

o o o 

coordinates of the observer, r , X and 0 are the position 

e e e 

coordinates of the emitting electron and R^, is the radius 
of the planet. 

e. Radiation from Electrons in a Plasma Medium 
The general formula for radiation in a plasma 
where there the observer is not in the medium is 


-T 


I = I (inc) e 
u> <u 


° + f 

Jo 


— T 

S (t) e dT 

U) 


where I is the radiation intensity at the frequency 


<JU 


0 ) , 1^ (inc) is the incident radiation upon the plasma medium. 
t is the optical dept and is the total optical depth. 

The quantity S^(t) under the integral sign is called the 
source function. 

For the quantity t, the optical depth, we have 


dT = -a ds 
0) 


( 26 ) 


(27) 


(28) 
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where a is the absorption coefficient and ds is the 

(D 


differential path length of the radiation. 


For the source function , we have 


^<d 2 a 

n u) 


r 

where n is the ray refractive index and J is the emission 
i* id 

coefficient . 

For J , the emission coefficient, we have 
w 


( 29 ) 


J = fj\ (P,a,\ , L) f (P ,a , X , L) 2rr sina p 2 d“ dP 

where P is the momentum of the electron and f(P,a,X o> L) 

is the distribution function of the electrons as a function 

of their momentum as well as a,\ ,L. However for a fixed 

o 

number of electrons N, we have 

2 

dN = f(P,a,X Q ,L) 2tt sina P da dP 
= p (E,X Q ,a ,L) dE da 

where E is the kinetic energy of the electron. 

The absorption coefficient a , is defined below as 


(30) 


(31) 


«. - - // V a ’V L ' E) TT 1 ^ slna p2 ® da (32) 

n ui a P 

r o 

For the ray refractive index n^, we have 
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2 

n = 
r 


2 • Q 

n sinG 


L 03 J 


2 n 1/2 


8 , c ° se + slne > 

r, ,Iin/l 1/2 

L 1 - + ( n ae 1 J 

03 


(33) 


In equation (33) r\ is the refractive index and 8 is the 
angle between the emission direction and the direction 
of the magnetic field. 

For the value of the refractive index, r\ , we would 

/V 

solve for the real part of the dielectric constant, u of 

/■N/ 

the medium. For u we have 

u = (n + i x) 2 (34) 

where x is the attenuation index. To solve for r\ , we 

« 

A/ 

would consider the Appelton equation for u. (see Heald and 
Wharton [ 13 ] . The Appelton equation is 


u = 


1 - 


2 2 
(«' /« ) 
E_ 


t 


2 2 2 
(<u /u) ) sin 0 
n . v b 

2,2 ,v 


2(1 - » p/ . - i - ^ 


+ 


4 4 

(U),/U3j sin 8 
b 4 


2.2 v 2 - 2 

4. (1.-03 /U) -1 ~) 03 

P 03 


03 

b 2 
+ — cos 6 


1/2 


1/2 


(35) 


where the + sign is for left circularly polarized waves 
and the - sign is for right circularly polarized waves. 


Also, oj is the observing frequency, o> b is the gyrofrequency 
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and v is the collision frequency. The plasma frequency 
is defined as 


di = 

P 


' o e ' 


1/2 


where n is the local thermal electron density of the 
e 

medium. 


( 36 ) 


f. Intensity vs Frequency for Various Latitudes 

Using Experimental Data on High Energy Electron 
Distributions in the Lower Van Allen Belt. 

We shall consider the intensity, I as a function 

of frequency, w for various observer latitudes at a height 

of 2 Earth radii. 

We shall require that the observer look along 
a line connecting the geomagnetic dipole center and the 
observer. 


For I , we have 

CD ’ 


S 


I 


U) 


(x e ,“) 


n 2 (S ) 
o 



J u> (S > U) > X e ) 
n 2 (S) 


dS 


l 

where n(S ) is the refractive index of the medium at the 
o 

observer and n(S) is the refractive index of the medium 
at the radiating source. I f we neglect the effects of 
the magnetic field, we get for n(S), 


n 2 (S) 


1. 


2 . 2 

tu /u) 

P 


(37) 


(38) 
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where o)^ is the plasma frequency of the medium at the 

' 2 

positions. For ou we have that 

P 

2 3 

«> = 3.1762 x 10 x N 

P t 

where is the thermal electron density at 'the position 
S. The values of N used for various values of S for 
these calculations were taken from data published by 
Vesecky [ 10 ] . 

The positions along r , the observer's distance 

o 

from the dipole center is determined by the formula 


S = L R_ cos X 
£ e 


where L is the L-shell containing the source and X is 

e 

the latitude of the source, 

The value of the emission coefficient J (S.m.x ) 

a) ’ ’ e 


is determined using the formula 

•n w (E,S,a,u),X ) p(E,S,a,X ) dE da 


V S -“- ,l e ) ' 


a 2 e 2 

'// 


a l E 1 


For the distribution function, p we will use the 
formula 

3 2 

cos x sin a 

p(E,S,a,X ) - 2— — N(S,E) 

(1+3 sin X > 1/4 
e 

where the position, energy distribution function N(S,E) is 
determined by using data compiled and published by Vette, 


( 39 ) 


(40) 


(41) 


( 42 ) 
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because the self-absorption of the electrons, determined 
by a^, has been showed to be very small by Vesecky [10], 
Considering figure (3) , we see that the intensity 
decreases as we increase the observer's latitude, This 

is to be expected from the nature of our electron distribution 
in terms of latitude X . Also, as we move up in latitude 
the upper limit on our region of radiation, L = 2. moves 
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deeper into the ambient plasma where the thermal electron 
density is greater. The mirroring of the electrons 
along field lines also contribute to this fall off in 
radiation for an increase in latitude. Notice that as 
we get up to the higher latitudes 35° and above the radiation 
falls off more rapidly than for lower latitudes as we 
pass the peak frequency. This is due to the fact that 
near these latitudes, the radiating medium between L = 1.35 
and L = 2. enters the atmosphere near the Earth’s surface. 

Now, all electrons near these latitudes must have pitch 
angles around 90 degrees or become lost in the lower 
atmosphere. This restriction causes the delta function 
S(Y) in the. formula for the coefficient of spontaneous 
emission equal to 

:6(r) = S (mu) - w (1 - 0 cos a cos 9)) (45) 

o 

to reduce to 

8(Y) = S(nw> - U)> , (46) 

o 

because a = 90 degrees. 

Considering equation (45) and (57), we see that as 
we increase the frequency u) for the low latitude case 
and the high latitude case, mio o in equation (46) must 
increase by a larger amount than for the low latitude case, 
equation (45) . This larger increase in harmonic number 
causes a corresponding larger decrease in intensity as 
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we look at higher frequencies. 


V. OBSERVED RADIO EMISSION 

We shall consider radiation received by the lower 
V antenna of the RAE-1 satellite. This is the antenna 
that is the primary receiver of the radiation coming from 
the direction of the lower Van Allen belt just below 
the satellite's orbit. 

Since the Van Allen belt is just below the satellites 

orbit we will assume it to be an extended source. Thus 

we shall assume that the intensity of the radiation 

received by the lower V antenna does not vary over the 

main antenna lobe for a given observer frequency i«, and 

satellite latitude X . This will allow us to replace 

o 

the brightness temperature, T fe by the antenna temperature, 

T . This gives the intensity, I (X ,u>) as 
A u) o 

- -h K % (47) 

8tt C 

where K is the Boltzman constant and c is the speed of 
light . 

Since the satellite is measuring the antenna tempera- 
ture at different latitudes and observing frequencies, we 
have that 

t a - t a <“.V (48) 
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In making comparisons between the radiation received 

by the RAE-1 lower V antenna and the calculated value 

of the radiation due to synchrotron radiation from our 

model, we shall multiply I (X ,u>) determined from (47), 

( 1 ) o 

by a factor of two. This is done because a linear type 
antenna only receives one degree of polarization of the 
radiation. 

We shall only consider the observing frequencies of 
1300 kHz and 2200 kHz for reasons stated below. Above 
2200 kHz radiation from the Earth's surface is able to 
penetrate the plasma medium of the lower magnetosphere 
and distort the quiet continuum emission. Below 1300 kHz 
the radiation is unable to penetrate the medium above the 
region where the majority of the radiation is being pro- 
duced, (at about L = 1.35). Thus, the majority of the 
radiation received below 1300 kHz is reflected radiation 
from the cosmic background. See figure (3). 

Figures (4) and (5) show a comparison of the plot 
of intensity, 1^ vs latitude, X (representing the average 
of twenty-five orbits of data) using the Ryle-Vonberg 
receiver data and a curve of intensity, I vs latitude 
predicted by the theoretical model for frequencies 1300 kHz 


18 



and 2200 kHz. Examination of these figures not. only 
show an increase in intensity as we approach the higher 
latitudes for the RAE lower V data curve in contrast to 
the theoretical curve which shows a decrease in intensity 
as we approach the higher latitudes but also show that 
the intensity at any latitude is much higher than that 
predicted by the model. The fact that the lowest intensity 
of the cosmic background observed by the RAE is greater 
than the peak intensity predicted by the model by a factor 
of about five made it virtually impossible for the lower 
Van Allen radiation belt emission to be detected using 
the present available data. However, if the present 
population of energetic electrons was increased to about 
five times its present value, the synchrotron emission 
might be observed. The above statement was made because 
it has been observed that during intense solar activity 
the population of the lower belt of energetic electrons 
is increased substantially. 

It has been showed experimentally that most of the 
electrons entering the polar regions of the magnetosphere 
have energies around 40 KeV [3]. We could not expect to 
obtain the quantity of synchrotron radiation observed by 
the RAE-1 in the auroral region (60 degrees latitude to 
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70 degrees latitude) from electrons with energy this 
low considering the cyclotron emission process. But if 
we assume that the velocity of the low energy electrons 
along the magnetic field lines is greater than the phase 
velocity of electromagnetic radiation in the plasma medium 
of thermal particles then we can consider the Cerenkov 
emission process. The energy transmitted by an electron 
to the plasma medium is 

u) E 

dE 1 p 2 r o -I 

as - 4 E “e 6 log [1 + K 1 


o 

where dE/dS is the energy per unit path length given up 

by the electron of energy E q to the excitation of plasma 

oscillations, m and e are the mass and charge of the 
e 

electron respectively, T is the temperature of the plasma 

medium and u) is the plasma frequency. Examination of the 
P 

above equation shows that the amount of energy transmitted 
is increased as the energy of the electron is decreased. 
The intensity, 1^ for electromagnetic Cerenkov radiation 
of low energy electrons in a magnetic field is 


(49) 


I = —7 HO " 31 u)J(P) l 
u) 2 tt 


1 Y - 1 ~ | watts 

0 2 (Y + X - l>J m 2. rad/sec . sr 


(50) 


where 


X 
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J(P) is the flux of the electrons and -t is the path 

length [9], Here u> is the plasma frequency, w is the 

P 

gyrof requency of the electron and <u is the observing 
frequency. 

For a large beam of particles the equation above for 
I no longer adequately describe the physical situation. 
Where as the first particle of the beam will be in a 
homogeneous plasma, the trailing particles will be affected 
by the plasma wave fields of the particles that precede 
them. These fields will move with the beam and so appear 
constant to the trailing particles within the Cerenkov 
cone. In this situation, the beam will become unstable 
and the particles will bunch. These bunched particles will 
radiate with partial coherence, greatly enhancing the total 
radiated energy. 

The condition for Cerenkov emission is 


— = cos 6^1 
nP 


(51) 


where P = .374 for 40 KeV electrons. For the refractive 
index, n we have 


1/2 


n = 11. + 


JL 


f(f b + f) 


where f , f, and f are the plasma, gyro and observing 
P b 

frequencies respectively in Hertzs. We must use the minus 
sign to insure n to be greater than unity for frequencies 


( 52 ) 
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I 

f not greater than the gyrofrequency , f fe . 

Let us choose the observing frequency f = 1.3 MHz 

and also choose the gyrofrequency very close to f. This 

will insure a large enough value of n to satisfy the 

Cerenkov condition for emission. 

For a latitude of 65°, the position where f is equal 

3 

to 1.3 MHz is at L-shelsJ. 1.2 or 2 x 10 Km above the Earth’s 

surface. The plasma frequency at this position is 284 kHz 

well below the observing frequency. 

Since f is much less than f, we have that 
P 

L z 1 * 1°' 37 f MHz J «> *■ [!• - V e 2 l < 53 > 

Now the observed electron flux density J(E > 40 Kev) is 
112 

4 x 10 electron/m .sec.sr [10], This gives for the 
intensity, 1^ 

~ . „-20! watts 

I = iO 5 

u) 2 

meter .sr.Hz 

which is of the order of magnitude of the radiation ob- 
served by the RAE-1 satellite. 

VI. CONCLUSION 

The model for gyro-synchrotron emission that we have 
developed will allow the spectral and spatial distribution 
of gyro-synchrotron emission from a celestial body containing 
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a dipole magnetic field to be computed given the distribu- 
tion of the trapped particles and the strength of the dipole 
field. 

We used an experimentally determined electron distribu- 
tion to compute what the spectral and spatial distribution 
of gyro-synchrotron emission should be for an observer 
viewing the inner radiating region of the Earth's lower 
magnetosphere (the lower Van Allen belt) at various 
positions . 

Comparisons were made with observations made with 
the RAE-1 satellite launched July 4, 1968. The satellite 
data showed that the cosmic background radiation is too 
large to allow for a direct observation of the Earth’s 
lower Van Allen radiation belt. 

It was also observed, upon examination of the RAE-1 
data that the emission from the lower magnetosphere tended 
to increase as the observing latitude was increased which 
is opposite to what our model for gyro-synchrotron emission 
predicts. However, in Section V, a calculation was made 
to obtain an estimate of what the emission should be from 
the 40 KeV electrons found in the higher latitudes considering 
the Cerenkov emission process. This calculation showed that 
the emission obtained from this process produced an intensity 
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of radiation of the same order of magnitude as was ob- 
served with the RAE-1 satellite. A more extensive inves- 
tigation of the Cerenkov effect as a possible mechanism 
responsible for the observed intensity from the lower 
magnetosphere might be fruitful. 

If one is allowed to speculate by saying that if, 
upon extensive investigation of the cosmic background 
emission received by the lower V antenna of the RAE-1 
satellite, a region of the sky is found where the cosmic 
background emission is lower than or nearly equal to the 
synchrotron emission predicted by the theoretical model 
we have developed, then it might be possible to detect the 
lower Van Allen radiation belt. 

Some specific examples of how the model developed 
in this paper might be useful are: 

1) the investigation of the magnetospheres of other 
planets (such as Jupiter which has a strong 
magnetic field) . Such an opportunity will present 
itself when NASA's Grand Tour Mission takes place 
and fly-bys of the outer planets of the solar system 
are made. 

2) The investigation of type IV sdlar bursts which 
are suspected to be caused by electrons spiraling 
in localized magnetic loops eminating from sun 
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spot regions of the sun's surface [7]. 

3) The investigation of cosmic background which 
is due to gyro-synchrotron emission caused 

by charged particles spiraling about magnetic 
field lines in inner steller space. 

4) The investigation of extra-galactic sources emitting 
synchrotron radiation, such as the strong radio- 
galaxies. 

The model developed by Vesecky [10] for synchrotron 
emission from the Earth's magnetic field confined the ob- 
server to the Earth's dipole equatorial plane, restricted 
the pitch angles of the electrons to 90 degrees and only 
allows for the computation of emission from paths confined 
to the dipole equatorial plane. 

The model developed by Chang [2] for synchrotron 
emission from the planet Jupiter does not allow for close 
up observation of the planet. It considers the planet 
to be a flat disk due to its great distance from the 
observer on the Earth. 

The model developed in this paper is more general 
than the two models stated above. It does not contain 


25 



their mentioned restrictions. However, it is restricted 
to the treatment of trapped particles in dipoie magnetic 
fields only. 


26 



BIBLIOGRAPHY 


[1] Bekefi, G. , "Radiation Processes in Plasmas", John 

Wiley and Sons, Inc., New York, 1966. 

[2] Chang, D.B., "Synchrotron Radiation from the Planet 

Jupiter", Boeing Scientific Research Laboratories 
DI -82-0060. 

’ [3] Davis, L.R., E. Berg, and L.H. Meredith, "Direct Measure- 

ments of Particle Fluxes in and near Auroras", Space 
Res., I, North-Holland Publishing Co., 1960, p. 721. 

[4] Katz, Ludwig, "Radiation Trapped in the Earth’s Magnetic 

Field", Edited by Billy M. McCormac, Gordon and 
Breach Science Publishers, New York (1966) 

[5] Razin, V.A., "Contribution to the Theory of Spectra of 

Radio Emission from Discrete Sources below 30 mc/s", 
News of Higher Educational Institutions, Ministry 
of Higher Education (USSR), Radio Physics Series 
3, 1960, pp. 73-97. 

[6] Stone, Robert G. , "Research Results from the Radio 

Astronomy Explorer", AIAA 6th Annual Meeting and 
Technical Display, AIAA, New York (1969) . 

[7] Stone, Robert G. , Private Communication (1970). 

[8] Van Allen, J.A. and L . A. Frank, "Radiation around the 

Earth to a Radial Distance of 107,400 km", Nature, 
183,430 (1959). 

[9] Van Allen, J.A., "Dynamics, Composition and Origin of 

the Geomagnetically-Trapped Corpuscular Radiation", 
Trans. IAU IX 99-136, (1962). 

[10] Vesecky, John F. , "Radio Frequency Synchrotron Radiation 

from Electrons Trapped in the Earth’s Magnetic 
Field”, SU-Sel-67-041 , 1967. 


27 



[llU Vette, J.I., A.B. Lucero, and J.A. Wright, "Models 

of the Trapped Radiation Environment, Volume II:. 
Inner and 0*ater Zone Electrons", NASA-SP-3024, 1966. 

[12] Weber, R.R., J.K. Alexander and R.G. Stone, "The Radio 

Astronomy Explorer Satellite" Goddard Space 
Flight Center, X-693-71-64. 

[13] Wharton, C.B. and M. A. Heald, "Plasma Diagnostics 

with Microwaves", John Wiley & Sons Inc. 

New York, 1965. 


28 





, INTENSITY( WATTS/M -sr Hz) 



4 5 6 7 8 9 

f (MHz) 

Graphs showing the intensity 1^ vs. frequency 
f, at various latitudes \ 0 , and a graph 
showing the total galactic background as a 
function of freauencv f 6 1 . 




(Watts/M 2 • sr • Hz) 



Figure 4 - A graph of I w (w , x 0 > vs \ Q for the case where 

experimental data is used to compute the expected 
radiation due to the synchrotron process compared 
with a graph of Io> (^ , X Q ) vs \ Q from the RAE-1 
lower V results. Here f = 1.3 MHz. 
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Figure 5 - A graph of Iuj(w,X 0 ) vs X Q for the case where 

experimental data is used to compute the expected 
radiation due to the synchrotron process compared 
with a graph of (u) , X Q ) vs X Q from the RAE-1 

lower V results. Here f = 2.2 MHz. 



